Strongly Reshaped Organic-Metal Interfaces: Tetracyanoethylene on Cu(lOO) 
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The interaction of the strong electron-acceptor tetracyanoethylene (TCNE) with the Cu(lOO) 
surface has been studied with scanning tunneling microscopy experiments and first-principles density 
functional theory calculations. We compare two different adsorption models with the experimental 
results and show that the molecular self-assembly is caused by a strong structural modification of 
the Cu(lOO) surface rather than the formation of a coordination network by diffusing Cu adatoms. 
Surface atoms become highly buckled and the chemisorption of TCNE is accompanied by a partial 
charge-transfer. 

PACS numbers; 68.43.Hn, 73.20.At, 68.43.Bc, 68.37.Ef 



Molecules containing multiple cyano (CN) groups 
such as tetracyanoethylene (TCNE) 7,7,8,8-tetra- 
cyanoquinodimethane (TCNQ) [2| or 2,3,5,6-tetrafluoro- 
7,7,8,8 tetracyanoquinodimethane (F4-TCNQ) 0,|3] rep- 
resent the archetype of strong electron acceptors. Fur- 
thermore, there is a promising development in the design 
of molecule-based magnets with high Curie temperatures 
that are based on strongly charged ligands anchored to 
paramagnetic transition metal atoms S S S| ■ How- 
ever, a significant amount of structural characterization 
still needs to be carried out in order to understand and 
describe the interaction of such strong acceptor ligands 
with metallic substrates. Recent studies on the adsorp- 
tion of TCNE and F4-TCNQ on single crystal surfaces 
such as Cu(lll) 0,13 and Au(lll) [ill^ reveal the exis- 
tence of strongly bound species where a large variation in 
charge density is invariably observed near the interface. 
In a comparative study of TCNE adsorption on various 
metal surfaces it was found that TCNE strongly self- 
assembles into chains on the Cu(lOO) surface. Moreover, 
scanning tunneling microscopy (STM) images contain ad- 
ditional bright protrusions and dark trench features that 
have been interpreted as consequences of a strongly re- 
constructed surface with buckled surface Cu atoms [lo| . 

In this Letter, we present results of first-principles den- 
sity functional theory (DFT) calculations and STM simu- 
lations as well as STM experiments of the structural and 
electronic properties of adsorbed TCNE on Cu(lOO). We 
found that TCNE is strongly adsorbed on the surface, 
and the adsorption gives rise to a reconstruction of the 
surface where the topmost Cu atoms bonded to TCNE 
are highly buckled with respect to the remaining surface 
Cu atoms. We have identified the presence of resonance 
peaks below the Fermi level that are associated to the 
highest occupied molecular orbital (HOMO) and to the 
lowest unoccupied molecular orbital (LUMO) of TCNE. 

Electronic structure calculations were carried out us- 
ing DFT within a local density approximation (LDA) [llj 
included in the siesta package For core electrons. 



norm-conserving pseudopotentials following a TrouUier 
and Martins scheme [l3| were used with relativistic cor- 
rections for Cu atoms. For the representation of va- 
lence electrons, we used an extended atomic basis set 
of polarized double-C type. The one-dimensional TCNE 
chain was built periodically from a supcrcell containing 
one molecule and a Cu(lOO) slab of three atomic lay- 
ers with about 3x8 Cu atoms in each layer. All atoms 
were allowed to relax except for Cu atoms in the bottom 
layer. Geometries were fully optimized using a Broyden 
scheme until a tolerance of 10""^ Ry/Bohr was reached 
for each system under study along with their deforma- 
tion and adsorption energy, total density of states (DOS), 
projected DOS (PDOS), MuUiken population analysis, 
and charge density. STM simulations were performed 
with our SPAGS-STM (Strongly Parallel Adaptive Grid 
Solvers - STM) software to evaluate topographic mode 
images and scanning tunneling spectra (STS). The soft- 
ware includes several algorithmic strategies such as par- 
allel computation of the tunnel currents [l^l and adap- 
tive grids that minimize the probing sites needed to ob- 
tain a high resolution image 15[. In STM simulations, 
the tunnel current was computed within a scattering ap- 
proach based on the Landauer-Biittiker formalism [l6| 
alon g w ith an extended Hiickel theory (EHT) Hamilto- 



man pj,ll7|. The sample preparation of TCNE/Cu(100) 
and STM experiments (using a Ft /It tip) were performed 
as described in Ref. Jjj. 

As a starting point, we intuitively placed the TCNE 
molecule in a highly symmetric adsorption site on a flat 
Cu(lOO) surface. As shown in Fig. [2a), the DFT-LDA 
geometry optimization reveals that a TCNE chain as- 
sembly induces an strong reconstruction of the topmost 
Cu(lOO) surface layer, most dramatically where the Cu- 
N bonds are formed. Within the one-dimensional chain 
structure, we can identify three different Cu atoms that 
are involved in the surface reshaping process: one Cu 
atom underneath the molecule chain is pushed below the 
surface plane by 0.3 A (marked as "1"), another surface 
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FIG. 1; DFT-LDA optimized structural geometry for the one- 
dimensional chain arrangement of TCNE on the Cu(lOO) sur- 
face (a) for the reconstructed model and (b) for the adatom 
model. A transparent plane is used to emphasize the magni- 
tude of reconstruction of the topmost Cu layer. 



atom (marked as "2") is highly buckled by 1.3 A out of 
the surface plane, and finally a third surface atom out- 
side the chain ( "3" ) opens a trench by sliding laterally by 
1.0 A toward the chain and 0.3 A out of the surface plane. 
The distance between highly buckled Cu atoms is 7.7 A 
along the chain direction and 6.6 A in the perpendicular 
direction. The energy difference between the buckled and 
a flat Cu(lOO) surface is 3.55 eV, which gives a deforma- 
tion energy of 1.78 eV per buckling. Strong interaction of 
adsorbates with substrates can induce nanoscale surface 
reshaping [1^]. For example, a single molecule CgoHgg 
at the step edge of a Cu(llO) surface was used to build 
a molecule-sized electrode [l9|. In another case, the re- 
moval of CeoHge self-assembly from a Cu(llO) surface 
left imprinted trenches underneath Metal surfaces 

roughening with Ceo were also reported [U, [g^, |2^ . 

On this reconstructed surface, TCNE is strongly 
bonded to Cu atoms. The energy needed to dissociate 
the TCNE/Cu(100) complex into individual but unrc- 
laxed components is 6.60 eV. Adsorbed TCNE is signifi- 
cantly deformed with respect to the gas phase structure, 
the relaxation of the molecule gives a deformation en- 
ergy of 0.71 eV. Hence, the formation of this strongly 
buckled structure would necessitate a minimal energy 
amount of (3.55 eV + 0.71 eV) = 4.26 eV, which is 
still 2.34 eV lower than the calculated dissociation en- 
ergy for the TCNE/Cu(100) complex. Although energies 
can be more accurately obtained by considering GGA 
functionals and basis set superposition error (BSSE) [^J] , 
the comparison of relative values obtained at LDA level 
should give an appropriate description. In addition, the 
deformation of isolated TCNE contributes to reduce its 
HOMO/LUMO gap from 2.61 eV to 1.98 eV. 

While the DFT-LDA optimization of the TCNE chain 
on Cu(lOO) suggests the described strong surface recon- 
struction, an entirely different self-assembly mechanism 
of molecules on Cu(lOO) is known to exist due to the fact 
that Cu adatoms from step edges can easily diffuse along 





FIG. 2: Topographic STM images of a TCNE chain (a) ob- 
tained ex per imentally with a Pt/Ir tip (V=l mV, 1=5 nA) 
from Ref. KLOl. and (b) calculated on a reconstructed Cu(lOO) 
surface with a Pt(lll) tip (V=l mV, 1=15 nA). (c) Calculated 
charge density of the valence electrons for the reconstructed 
model (a black plane is used to highlight the density drop), 
(d) Calculated STM image of a TCNE chain for the adatom 
model. All image domains are 31 A x 18 A. 



the Cu(lOO) terraces ^25*]. These Cu adatoms can connect 
to molecules and promote ordered molecular structures 
via formation of a coordination network or- 
der to consider such a self-assembly process for the TCNE 
chains on Cu(lOO), we also performed DFT-LDA calcula- 
tions on the adsorption of TCNE with four Cu adatoms 
located in four-fold sites of the Cu(lOO) surface. The cal- 
culated optimized geometry is shown in Fig.[ljb). In this 
case, the surface layer remains unreconstructed, and the 
adatoms have left their four-fold positions to establish co- 
valent bonding with CN groups. The surface structure of 
this adatom model shows a distance between Cu adatoms 
of 7.7 A along and 5.4 A across the chain, the latter being 
significantly shorter (by 1.2 A) than in the reconstructed 
model. In addition, the dissociation energy for this com- 
plex is 7.67 eV, which makes the TCNE more tightly 
bound to the surface than in the reconstructed model. 

Although the exclusion of one of the two mechanisms 
remains a difficult task without investigating the dynam- 
ics of the processes, we can try to discriminate between 
them by comparing how well simulated STM images of 
both models resemble the experimental STM images of 
the TCNE chains, which show three characteristic fea- 
tures (Fig. [21(a)): (i) The TCNE molecules appear as 
a bright, elongated protrusion centered around the C=C 
bond and four faint lobes close to the CN groups (referred 
to as short legs in Ref. 10). {ii) Additional protrusions 
next to CN groups of neighboring molecules above and 
below the TCNE chain {extended legs) are presumably 
caused either by buckled Cu surface atoms or by bound 
Cu adatoms [30|. {in) Dark trenches surround the buck- 
led atoms or adatoms above and below the TCNE chain. 

The STM simulation of the reconstructed model 
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probed with a Pt(lll) tip is displayed in Fig. HJb) 
The theoretical STM image can reproduce the first two 
features, although the short and extended legs appear 
more blurred. The third feature, the dark trench, is faint 
but visible in the STM simulation. Moreover, we observe 
a clear drop in the calculated charge density map close 
to the buckled atoms (Fig. HIc)), which can be associ- 
ated to the dark trench area. Overall, the characteristic 
experimental features are reproduced fairly well. In com- 
parison, the simulated STM image of the adatom model 
is shown in Fig.[2I^d). While also this model qualitatively 
reproduces the bright TCNE center as well as the short 
and extended leg features (that also appear blurred), no 
dark trench can be reproduced at all. Hence, a main 
characteristic feature of the experimental images is miss- 
ing in this model. 

Finally, we quantitatively compare the features re- 
lated to the extended legs in experimental and simulated 
STM images (see arrows in Fig. [2|). From the experi- 
ments we find distances of 7.7±0.1 A along the chain 
and 7.3±0.3 A across. Both theoretical models show a 
distance of 7.7±0.2 A along the chain (i.e., three times 
the Cu nearest- neighbor distance) , which is in very good 
agreement with the experiment. Across the TCNE chain, 
the situation changes: for the reconstructed model, the 
simulated STM image shows a distance of 7.7±0.2 A, 
which is slightly larger than the actual distance of the 
highly buckled atoms, because the Cu atom "3" next to 
the trench has a significant contribution to the apparent 
height in the STM image. The simulated STM image of 
the adatom model exhibits a distance of only 5.3±0.2 A 
across the chain, similar to the actual adatom distance 
across the TCNE chain. This is in strong disagreement to 
the experiment. Consequently, the reconstructed model 
reproduces the experimental observations far better, and 
thus we conclude that the strong buckling of the Cu(lOO) 
surface is most likely the driving mechanism for the ob- 
served self-assembly of TCNE molecules into chains. In 
addition, the DFT and STM calculations agree to sup- 
port that mobile adatoms do not apparently participate 
in the formation of the one-dimensional TCNE chain: 
the most stable geometry is not experimentally observed 
by STM and the STM simulations on the reconstructed 
model are accurately reproducing the experimental STM 
results. 

Another important aspect of the interaction of TCNE 
with surfaces is the magnitude of charge transfer. The 
formation of TCNE^^ or TCNE^~ species is generally as- 
sumed when low work function metals such as Na or Ca 
are used [ij . The situation is slightly different with tran- 
sition metals for which the accumulation of net charge on 
TCNE depends on the magnitude of TCNE^metal do- 
nation but also on the metal— >TCNE backdonation 
Hence, a fractional net charge on TCNE would reflect 
that a charge transferred to the metal (donation) is sig- 
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FIG. 3: (a) Total DOS and projected DOS for a buckled Cu 
surface without TCNE and (b) after adsorption of TCNE, 
and total DOS of an isolated TCNE one-dimensional chain, 
(c) Wavefunction contours of the molecular orbitals located 
at -8.9 eV and -5.5 eV in (b) and HOMO and LUMO of 
isolated TCNE. 



nificantly returned back to the molecule (backdonation) . 

Fig. [3ja) shows the total DOS and various PDOS for 
a buckled surface without the TCNE molecule. The to- 
tal DOS and PDOS over most Cu atoms of the buckled 
surface are similar to the DOS diagram usually observed 
for this Cu surface [1^, except for highly buckled Cu 
atoms for which the states are more strongly localized 
(sharp and intense) and more weakly bonded than other 
Cu surface states. The adsorption of TCNE is considered 
in Fig. ^h) : the sharp peak of highly buckled Cu atoms 
at —6.0 eV is now smeared out within the broad PDOS 
peaks at higher binding energy. The disappearance of 
the localized surface states into the higher binding en- 
ergy manifold suggests an electron depletion of buckled 
Cu atoms and a strong mixing of states with TCNE or- 
bitals. In Fig. [2Ib), the DOS of an isolated chain shows 
quite narrow peaks while the PDOS of adsorbed TCNE is 
much broader. Therefore, such broadening results from 
the important mixing between TCNE and Cu(lOO) states 
rather than from the coupling between TCNE molecules 
in the chain. 

The nature of two TCNE-related molecular orbitals 
near the Fermi level are shown in Fig. Efc). We can 
easily attribute the peaks at —8.9 eV and —5.5 eV of 
adsorbed TCNE to the HOMO (-8.1 eV) and LUMO 
(-6.1 eV) of isolated TCNE, respectively. The shifts 
with respect to the isolated TCNE chain can be ex- 
plained in terms of charge transfer: the HOMO moves to- 
ward high binding energies (electron donation to Cu) and 
the LUMO moves toward low binding energies (electron 
backdonation to TCNE). The contribution from TCNE 
to the wavefunction of the TCNE/Cu(100) complex cen- 
tered at —5.5 eV is ~90%, i.e. TCNE accepts at least 
1.8 |e| from Cu(lOO). MuUiken population analysis of the 
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TCNE/Cu(100) system gives a net charge of 0.30 |e| on 
TCNE, which indicates that Cu receives ~1.5 |e| from 
HOMO and HOMO-n states of TCNE. A comparable 
charge transfer of 0.55 |e| from a polycrystaUine Cu sur- 
face to TCNE has been deduced from surface-enhanced 
Raman scattering experiments [2^. The adsorption of 
the electron acceptor F4-TCNQ on the Cu(lll) surface 
also leaves a net charge on the molecule of approximately 
0.6 |e| i. 

In summary, slab DFT-LDA calculations, STM ex- 
periments and simulations reveal that the adsorption of 
TCNE molecules on Cu(lOO) is accompanied by a re- 
construction on the topmost Cu surface layer. By com- 
parison, we exclude the possibility of an alternative self- 
assembly mechanism based on diffusing Cu adatoms. Re- 
constructing Cu atoms that bind directly to the CN 
groups of TCNE are strongly buckled out of the surface 
plane by 1.3 A. The presence of dark trenches in STM ex- 
periments and simulations is attributed to a drop in the 
surface charge density that is caused by the concerted lat- 
eral displacement of Cu atoms contributing to the buck- 
ling process. TCNE interacts strongly with the Cu(lOO) 
surface leading to a large adsorption energy. However, 
the net charge transfer onto the TCNE moiety is rela- 
tively weak due to a significant amount of charge that is 
backdonated to the molecule. 
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